Optically-Induced Polarons in Bose-Einstein Condensates: Monitoring Composite 

Quasiparticle Decay 
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Nonresonant light-scattering off atomic Bose-Einstein condensates (BECs) is predicted to give rise 
to hitherto unexplored composite quasiparticles: unstable polarons, i.e., local "impurities" dressed 
by virtual phonons. Optical monitoring of their spontaneous decay can display either Zeno or anti- 
Zeno deviations from the Golden Rule, and thereby probe the temporal correlations of elementary 
excitations in BECs. 
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The evolution of a quantum state coupled to a con- 
tinuum is coherent and effectively reversible over the so- 



called "correlation time" t c 



This time marks the du- 



ration of the non-exponential initial stage of decay of the 
unstable state [1-3], or, equivalently, its response time to 
fast (impulsive) perturbations. 

Little is known about the short-time dynamics of 
impulsively-perturbed quantum many-body systems. 
Atomic Bose-Einstein condensates (BECs) are especially 
suitable for the exploration of such effects, since their 
elementary excitations, which are describable as quasi- 
particles [4-10], may remain correlated for a long time. 
Here we consider the short-time dynamics of a compos- 
ite quasiparticle, consisting of a moving "impurity" atom 
(an atom in a ground-state sublevel different from the 
rest of the BEC) "dressed" by a cloud of virtual phonons 
due to deformation of its vicinity. This quasiparticle may 
be viewed as the BEC analog of the solid-state polaron 
[11]. The time necessary for dressing an impurity atom 
by phonons is just the correlation time of the quasiparti- 
cle, which is shown to be exceedingly long (~ 1 ms). The 
polaronic effect in a BEC is a genuine example of the 
response of a many-body quantum system to the pres- 
ence of a probe particle. Direct measurement of t corr is 
a formidable task, since bare and dressed impurities are 
hardly distinguishable spectroscopically. 

We therefore propose an indirect method of measuring 
tcorr'- optically- induced formation of polarons by Raman 
scattering followed by real-time optical monitoring of the 
products after their decay. The decay products are pair- 
correlated elementary excitations that are unambiguously 
detectable. Their time-resolved monitoring can reveal 
tcorr, which may be manifest by either the Zeno or the 
anti-Zeno effects [1-3], namely, slowdown or speedup of 
the decay rate compared to its Golden Rule value. Fluc- 
tuations of the laser fields, on time scales comparable 
to the very long t corr of the BEC, in the ms range, can 
be used to interrogate these features more easily. Simi- 



lar results are obtained for the Bragg process, where, in 
contrast to the Raman process, the atomic internal state 
is not changed. 
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FIG. 1. (a) The Feynman diagram for the formation and 
decay of the polaron includes the Raman beams (curved 
lines), the unstable polaron (double dashed line) and the cor- 
related products: the bare impurity (dashed line) and BEC 
Bogoliubov excitation (solid line) (b) The momentum distri- 
bution, (parallel and perpendicular to Tiq, which is indicated 
by =>) of the impurity atoms (dark grey) and Bogoliubov 
excitations (light grey) of the BEC (central dark spot) pro- 
duced by Raman off-resonant scattering with hA — 0.66pi, 
hq = 0.14mc s . The momentum shells are broadened due to 
finite pulse duration (120h/fi). These shells should be ob- 
served in time-of-flight images of the excited condensate. 

We shall primarily consider polaron formation by Ra- 
man scattering off an atom in the BEC, as shown in Fig. 1 
(a). Two laser beams are arranged such that momen- 
tum transfer of hq is accompanied by energy transfer 
of hA to the atom. The effective Rabi frequency asso- 
ciated with the induced Raman two-photon transition is 
CI. Proper choice of the laser-beam frequencies and polar- 
izations ensures that the atom is transferred to another 
sublevel of the ground state, thereby leaving the con- 
densate and forming an impurity atom [12]. The energy 
of such an impurity atom is given by its kinetic energy, 
h 2 q 2 /(2m), plus the mean-field shift, 4arh an/m, a be- 
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ing the scattering length for a collision between the im- 
purity and condensate atoms, plus the energy-difference 
Ed between the two (Zeeman and/or hyperfme) sub- 
levels involved. We can incorporate the difference be- 
tween the mean-field shift and the chemical potential 
fi = Airh 2 aon/m, a being the BEC scattering length, as 
well as Ed into a new definition of the Raman-transferred 
energy, HA = HA + Airh 2 (cto — a)n/m — Ed- This enables 
us to identify the energy of the impurity atom with its 
kinetic energy. We shall assume large blue two-photon 
laser detuning, A > hq 2 / (2m), O < A, so as to avoid the 
need to deal with strongly-driven atoms [13]. 
The relevant Hamiltonian for a uniform BEC is 



H — H at + H int , 



(1) 



the first term describing the atomic system itself, and 
the second one — its interaction with the laser radiation. 
These terms are: 

h 2 k 2 



(e-ubt + gkfok) , (2) 



H mt = hilVN (e- lAt pl + e jAt /3 q ) . (3) 



Here the creation and annihilation operators of the 
BEC elementary Bogoliubov excitations are denoted by 
& k , &k, and those of impurity atoms by /3k, = 
ky/\hk/ (2m)] 2 + c 2 is the frequency of the BEC elemen- 
tary excitation with the momentum %k, c s — \J fi/m is 
the speed of sound in the BEC, Sk = hk 2 / (2muJk) is the 
BEC static structure factor, = /3k + k'/^ k ' i s * nc 
operator of the impurity momentum shift, and V is the 
quantization volume. 

Now we define the dressed (polaronic) states as eigen- 
states of the atomic-system Hamiltonian (2) and express 
the field-atom interaction Hamiltonian (3) in the basis 
of these states. We omit the intermediate calculations, 
which are quite involved but straightforward [11], and 
present the result: 



H, 



int 



n 

Vv 



■Mq,k(i)Cq- k ,kCO+H.C. 



(4) 



where cq is the operator of annihilation of a bosonic 
atom in the BEC state, and Cq_ k k is the operator of 
creation of a correlated pair consisting of a dressed im- 
purity atom with momentum h(q — k) and a elementary 
Bogoliubov excitation with momentum hk. The latter 
operator obeys the usual bosonic commutation rules as 



long as the population of the state |q— k, k)d = c 



t 

q-k,k 
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is much less than 1. The energy of the state |q — k, k)^ 
is e | q -k|,fc = + h 2 (c[ — k) 2 /(2m), provided that the 
energy correction (to second order in a), which has to 



be calculated using the scattering- length renormalization 
[5], is included in the definition of HA. The matrix ele- 
ment of the off-resonant two-photon transition that cou- 
ples the initial vacuum state |0) to the state |q — k, k)^ 
is 



-iAt 



m[ujk + hk 2 /(2m) — 7iqk/ra] 



(5) 



The relevant Feynman diagram is shown in Fig. 1 
(a). The operator Heisenberg equations obtainable from 
Eq. (4) are 



d 



c = ^-M q , k c q _k,k, 



d 



''77-Cq-k,k = e|q-k|,feC q -k,k + Mq y kC - 



(6) 



(7) 



Equations (6, 7) represent the well-known model of a 
discrete state coupled to a continuum. Assuming a per- 
turbative production fraction, the time-dependent rate 
of creation of correlated pairs, consisting of an impurity 
atom and a Bogoliubov excitation is 



r(t) 



d 3 k |2 2simVi) 



(27T) 



(8) 



where v = A — ujk — ft(q — k) 2 /(2m). 

In what follows we assume that hq -C mc s , i.e., 
the virtual impurity atoms move at a subsonic velocity. 
Whereas the detuning of the beams A from the Raman 
resonance, ensures a significantly larger momenta (and 
subsequent energy) for the decay products. 

The quantity ?ikq/m can then be neglected compared 
to u>k, and Eq. (8) can be approximately expressed as 



r(t) 



=r *r 



■ rin[(f-A,)t] 

tt(<5 - A g ) v ; 



Here 



4^7? 




(9) 



(10) 



is the characteristic decay rate, G(S) = (hS/fj,) [ (h5/n) + 
1 ]~ 5 / 2 is the dimensionless response function of the uni- 
form BEC, and A q = A - hq 2 /(2m) w A. The in- 
tegral in Eq. (9) can be performed analytically. The 
proportionality of the rate prefactor T* in Eq. (10) to 



V na 3 clearly indicates that the process under consider- 
ation relies on non-mean-field (impurity scattering) ef- 
fects. For fit/U 3> 1 Eq. (9) gives the Golden Rule rate 
T GR = r,G(A,). 

In Fig. 2 we display T(t)/ToR for three different de- 
tunings. For fit/h^A the drastic deviations from T GR 
imply that we are within the correlation time t corr of the 
elementary excitations before irreversibility sets in. 

The Raman excitation fraction, namely, the num- 
ber of correlated pairs created per BEC atom, P(t) = 



2 



Jg dt' r(i'), is assumed to be small so as to avoid conden- 
sate depletion: F< 1. Then Eq. (9) yields 



pa) 



dS 



2 sin 2 [((5- A q )t/2] 
tt(5 - A q ) 2 



G(8). (11) 




FIG. 2. The creation rate V (normalized to Tgr) of the 
laser-induced impurity atoms and Bogoliubov excitations as 
a function of dimensionless time for hA q //j, equal to 2.0 (solid 
line), 0.66 (long-dashed line), and 0.07 (short-dashed line). 
Typical correlation time is t corr . 
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FIG. 3. The creation rate r T (oo) of laser-induced impurity 
atoms and Bogoliubov excitations versus the dephasing time 
r for the same hA q /fi as in Fig. 2. AZE and QZE regions are 
indicated. Inset: the medium response function G(S) (thin 
line) and typical modulation spectrum F(6) (thick line), both 
in arbitrary units. 

The long correlation time of the BEC can lead to a 
significant deviation of P{t) from that expected from 
Golden Rule. In principle, this deviation is observable, 
even at long times, as a shift from the Golden Rule pre- 
diction, P(t) = Tent. However, this deviation is difficult 
to observe for realistic experimental parameters, as we 
show below. 

We therefore suggest an alternative scheme, which al- 
lows us to observe large deviations from the Golden Rule 
rate, even at the limit of large times. Suppose that the 
difference of the two Raman laser frequencies is randomly 
modulated [3], so that the well-defined detuning A q is 



replaced by the spectral distribution F(S) normalized 
to 1. Its spectral r.m.s. fluctuation (around the mean 
frequency A q ), t _1 , may be regarded as the inverse de- 
phasing time. Under these conditions, even the long-time 
limit of the excitation rate r T (oo) obtained at t ^> fij \x 
may strongly differ from Tor- As follows from the uni- 
versal formulas for the QZE [3], Eqs. (9, 11) take in this 
case the form 



r T (t) = r. 



p T (t) = r* 



d5 



fOO PC 

/ ds 

Jo Jo 



dS' 



dS' 



sm[(6-6')t] 

ir(S-S>) 
2 sin 2 {(5 -S')t/2] 
n(S-S') 2 



G(5)F(S'), (12) 



G(S)F(5'). 



Equation (12) implies that the modulation spectrum 
F(S) controls the asymptotic decay rate r r (oo) and ex- 
citation fraction P T (oo), so that a suitable choice of the 
modulation spectrum may yield conspicuous deviations 
from the Golden Rule, in contrast to the case of Eq. (11). 

Specifically, we choose the modulation spectrum F(S) 
with a sharp low-frequency cut-off in order to suppress 
undesired transitions from |0) to the state containing only 
a dressed impurity with the momentum hq and no real 
Bogoliubov excitations. These requirements are satisfied 
by F(S) = C5~ 3 / 2 for Si < 5 < S 2 and zero otherwise (see 
inset of Fig. 3). The coefficient C ensures normalization 
to 1 , and the cut-off frequencies Si , S2 are chosen so that 
the mean frequency is A q and the r.m.s. fluctuation is 
t _1 . The resulting values for r r (oo) are shown in Fig. 3 
as functions of r for different values of A q . 

The convolution of G(S) and F(S) may give rise to 
either the Quantum Zeno effect (QZE) or the opposite, 
anti-Zeno effect (AZE), i.e. decay speedup at short times 
[2]. As expected [1-3], the QZE always takes place for ex- 
tremely small times r. For large detunings, the modula- 
tion spreads the accessible final states towards the range 
where the response G(S) increases as r" 1 grows, and the 
irreversible transitions are accelerated (AZE takes place) 
[3] . Alternatively, for small and moderate detunings, the 
modulation spreads the final states over a range of ener- 
gies such that G decreases as t _1 grows, thereby bringing 
about the QZE. Comparison of Figs. 2 and 3 indicates 
that the dephasing time scale r associated with the QZE 
can be identified with the correlation time t corr of ele- 
mentary excitations. 

We obtain qualitatively similar results for off-resonant 
Bragg scattering, although there are some subtle caveats, 
due to the indistinguishability of the excitations, as op- 
posed to the Raman process. The general expression for 
the rate Tsraggit) of this process is quite cumbersome. 
Here we present the t — > 00 limit of V Bragg{oo) for the 
particular case A » /j,/K ^> hq 2 /(2m): 

n \ 3 / 2 hfi 2 
,ka) (j, 

Remarkably, r Bragg docs not depend on q, as long as 
q£> 1, I being the length of the trapped BEC sample. 



(13) 
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The Bragg process is characterized by a correlation time, 
which is similar to that of the Raman process and, hence, 
can be also revealed via the QZE and AZE. 




t [ms] 

FIG. 4. The Raman excitation fraction (dimensionless), as 
a function of time for a S7 Rb BEC with n = 4 • 10 14 cm~ 3 
and = A q « 1.3 • 10 4 s" 1 . The solid line is the Golden 
Rule result. The dashed line is the prediction for a cw Ra- 
man experiment, and the dotted line is the prediction for a 
frequency modulated Raman experiment, with the laser beam 
correlation time r = 1.0 h/fi. The modulation leads to a clear 
asymptotic deviation from the Golden Rule result that is eas- 
ily observed experimentally, in contrast to the case of a cw 
Raman excitation. 

We used for our plots the response function G(5) of 
a uniform BEC. It can be generalized to the local den- 
sity approach if the length of the trapped BEC sample 
satisfies q£ > 1 (usually £ ~ 10~ 3 cm). This inequal- 
ity is compatible with the subsonic limit for hq/m under 
typical experimental conditions [12,14]. 

As an example, we consider a 87 Rb BEC with n = 
4 • 10 14 cm" 3 [8] and Q = A q w 1.3 • 10 4 s"\ for the 
Raman process. We expect, after the excitation, to ob- 
serve by the time-of-flight technique [14] unambiguous 
momentum shells of Raman and phonon excitations, as 
shown in Fig. 1 (b). Although the probability of di- 
rect excitation of dressed impurities with momentum hq 
may be non-negligible because of broadening of various 
kinds (inhomogeneous, Fourier etc.), impurity atoms pro- 
duced via this channel are well-separated by energy and, 
therefore, distinguishable by a time-of-flight measure- 
ment from those produced by the off-resonant Raman 
process, as indicated by the double arrow in Fig. 1 (b). 

In Fig. 4 we show the total excitation fraction P(t) 
(dashed line), which should be contrasted with the 
Golden Rule production of excitations (solid line), in- 
dicating the presence of the QZE. Experimentally, mea- 
suring the population in the excitation shells may allow 
us to resolve t corr by observing a shift in the asymptotic 
increase of excitations. However, this is quite challeng- 
ing. Therefore, we propose to use the frequency mod- 
ulation scheme described above. This will lead to an 
easily measured deviation of the excitation fraction from 
its Golden Rule counterpart for all times, as shown in 
the dotted line, calculated for a laser beam correlation 
time r = 1.0 ti/fi. Specifically, after 1 ms this deviation 



is nearly 40%, and the corresponding excitation fraction 
P T i=a 0.05 is easily detectable experimentally [8,12]. 

To summarize, we have presented a theory of off- 
resonant stimulated Raman and Bragg scattering in 
BECs, followed by the hitherto unexplored creation of 
correlated pairs of matter-wave quanta - the decay of 
a composite polaron. These processes are found to be 
highly suitable for experimental observation of the quan- 
tum Zeno and anti-Zeno effects. These effects may be 
used as unique probes of correlation times of the BEC 
response to fast (impulsive) perturbations, and the onset 
of irreversibility in the quasiparticle formation. 
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